Magnetotransport and thermoelectricity in disordered graphene 
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( ^^ ) We have studied the electric and thermal response of two-dimensional Dirac-fermions in a quan- 

tizing magnetic field in the presence of localized disorder. The electric and heat current operators 
in the presence of magnetic field are derived. The self-energy due to impurities is calculated self- 
C"| , consistently, and depends strongly on the frequency and field strength, resulting in asymmetric 

■ peaks in the density of states at the Landau level energies, and small islands connecting them. 
The Shubnikov-de Haas oscillations remain periodic in 1/B, in spite of the distinct quantization of 

, quasiparticle orbits compared to normal metals. The Seebeck coefficient depends strongly on the 

■ field strength and orientation. For finite field and chemical potential, the Wiedemann-Franz law 
can be violated. 

^ | PACS numbers: 81.05.Uw,71.10.-w,73.43.Qt 
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C« . I. INTRODUCTION 

& ■ 

Recent advances of nanotcchnology have made the creation and investigation of two dimensional carbon, called 
graphene, possible^ 2 -! 3 -^. It is a monolayer of carbon atoms packed densely in a honeycomb structure. In spite of 
*^ , being few atom thick, these systems were found to be stable and ready for exploration. One of the most intriguing 
£h ' property of graphene is, that its charge carriers are well described by the relativistic Dirac's equation, and are 
O , two-dimensional Dirac fermions^. This opens the possibility of investigating "relativistic" phenomena at a speed of 
i ~ 10 6 m/s (the Fermi velocity of graphene), l/300th the speed of light. The linear, Dirac-like spectrum causes the 

\ density of states to increase linearly with energy, which is to be contrasted the constant density of states of normal 

■ metals. Due to this peculiar property, the response of graphene to external probes is expected to be unusual. This 
manifests itself in the anomalous integer quantum Hall effect, which occurs at half-integer filling factors, and in 
the presence of universal minimal value of the conductance. The dependence of the thermal conductivity on applied 
magnetic field has been measured in highly oriented pyrolytic graphite 7 -^. 

Dirac-fermions show up in other systems, at least from the theoretical side. They characterize the low energy 
properties of orbital antiferromagnets, a density wave system with a gap of d-wave symmetry 9 -^. A similar model has 
f***- | been proposed for the pseudogap phase of high T c cuprate superconductors, known as d-density wave, with peculiar 
electronic properties^. A similar system was also mentioned in the context of heavy fermion material URu2Si2, which 
' shows a clear phase transition at 17 K without any obvious long range order, detectable by X-ray or NMR experiments. 
Its low temperature phase was attributed to another spin density wave with a d-wave gap 12 ' 13 Experimentally, the 
aforementioned materials possess unusual electric and thermal responses as a function of temperature and magnetic 
^ i fields. 

Therefore, the interest in studying the transport properties of two-dimensional Dirac-fermions is not surprising. 
Q Gusynin, Sharapov and coworkers have studied exhaustiveb y 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 the electric and thermal response of two- 
O ■ dimensional systems with linear energy spectrum, with special emphasis on the Wiedemann-Franz law and magnetic 
oscillations. However, their self energy due to scattering from impurities was not determined in a self-consistent 
manner, but rather they assumed a constant, energy, magnetic field and temperature independent scattering rate. 
Moreover, they completely neglected the real part of the self energy, responsible for the shift of energy levels. Nev- 
ertheless, they derived beautiful analytical formulas for the various transport coefficients, which, although suffering 
from the above limitations, turned out to be useful in explaining experiments 6 -. 

Impurity scattering can be taken into account in the presence of quantizing magnetic field in the usual self-consistent 
way^ 3 .. This program has been carried out, among many others 2 ^, by Peres et alj2£. In their work, the full self- 
consistent Born approximation was used before taking the strength of the impurity potential to infinity. They studied 
the frequency dependence of the electric conductivity for various fields, but never entered into the realm of thermal 
transport. Parallel studies have also been performed in the limit of wcak-scatterer o 26 ' 27 . 

In this paper, we extend the work of Refs. [l6llT7lfl8lfl9l and determine self-consistently the energy and magnetic 
field dependent self-energies and study the Seebeck coefficient as well, and also generalize Ref. [25] to include thermo- 
electricity. We study Dirac-fermions in a Landau quantizing magnetic field (B) in the presence of scatterers, allowing 
for arbitrary field orientations. In a way, our study here bridges between the efforts of the previous groups. After 
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the introduction of the general formalism, we determine the electric and heat current operators, essential for further 
steps. By introducing impurities in the system, we can study the quasiparticle density of states, the electric and 
heat conductivity, the Seebeck coefficient and the Wicdemann-Franz law as a function of magnetic field strength and 
orientation and temperature. For high fields, the discrete nature of the Landau levels is revealed in the density of 
states in the form of asymmetric peaks at Landau level energies (far from being Lorentzians), which smoothen with 
decreasing field. Shubnikov-de Haas oscillation are visible in all transport coefficients, periodic in 1/B, similarly to 
normal metals^. The angular dependent conductivity oscillations become more pronounced with increasing field. 
The chemical potential dependence of the conductivity resembles closely the experimental findings^. The Seebeck 
coefficient depends strongly on the applied field and temperature. 



II. LANDAU QUANTIZATION, ELECTRIC AND HEAT CURRENT 

The Hamiltonian of non-interacting quasiparticles living on a single graphene sheet is given b y 25 ! 29 ' 30 : 

H = -v F ^ Hfy + eA i(r)) > (!) 

where a/s are the Pauli matrices, and stand for Bloch states residing on the two different sublattices of the bipartite 
hexagonal lattice of graphen o 18 ' 25 . Strictly speaking, the Hamiltonian above describes quasiparticles around the K 
points of the Brillouin zone, where the spectrum vanishes. The vector potential for a constant, arbitrarily oriented 
magnetic field reads as A(r) = (— By cosd, 0, B(y sva.9 cos <f> — x sin 9 sin (f>) ) , where 9 is the angle the magnetic field 
makes from the z axis, and <p is the in-plane polar-angle measured from the a;-axis. We have dropped the Zeeman 
term, its energy would be negligible with respect to energy of the Landau levels, Eq. using vp ~ 10 6 m/s, 

characteristic to graphene. Eq. ([T]) applies for both spin directions. 

In the absence of magnetic field, the energy spectrum of the system is given by 

£7(k)=±«*.|k|. (2) 

This describes massless relativistic fermions with spectrum consisting of two cones, touching each other at the end- 
points. From this, the density of states per spin follows as 

pH = -V^-B(k)) = ^ / 'kdk6(u±v F k) = ^, (3) 
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where k c is the cutoff, D = vpk c is the bandwidth, and A c = 47r/fc 2 is the area of the hexagonal unit cell. 

In the presence of magnetic field, the eigenvalue problem of this Hamiltonian (Hq^> = Ety) can readily be solved 2 ^. 
For the zero energy mode (E=0), the eigenfunction is obtained as 



and the two components of the spinor describe the two bands. The energy of the other modes read as 



E(n, a) — au> c \/ n + 1 (5) 



with a = ±1, u> c = VFy/2e\B cos(6)\ is the cyclotron frequency, n = 0, 1, 2, Only the perpendicular component 

of the field enters into these expressions, and by tilting the field away from the perpendicular direction corresponds 
to a smaller effective field. The sum over integer n's is cut off at N given by N + 1 = (D/ui c ) 2 , which means that 
we consider 2N + 3 Landau levels altogether. For later convenience, we define a magnetic field _B , whose cyclotron 
frequency is equal to the bandwidth (u> c = D). 
The corresponding wave function is 

* , (v)- e —( Mv-kl%) \ 

with cyclotron length If, = lyeB. Here <fi n (x) is the nth eigenfunction of the usual one-dimensional harmonic oscillator. 
The electron-field operator can be built up from these functions as 



*(r)=E 



k L 



(7) 



FIG. 1: The structure of the Landau levels is visualized schematically for the first few levels. 



The Green's functions of these new operators do not depend on k, and read as 

G (iu! n , k) = - — , 

1 

G {iui n ,k,n, a) 



iu n - E(n,a) 



(8) 
(9) 



for Cfe and Ck, n .a, respectively, and u> n is the fermionic Matsubara frequency. 

With the use of these, we can determine the electric and heat current operators of the system. Following Mahan^L, 
we define the polarization operator as 



¥ = \J dr M r ) + <°( r ) r ) 



(10) 



with p(r) = ^E' + (r)^E'(r) giving the charge density, and the symmetric combination ensures hermiticity. The total 
current is its time derivative, which follows as 



3 = dtP = i[H,-p}. 

By performing the necessary steps, after straightforward calculations this yields^ 



J x = v F e^2 



p , a 



J y = ivpe 



P ■ <x 



~7k\ C P C P,0,a + C p,0..a C PJ + 2-^ 9 '^n+l.a^"^ + Cp, n ,\ C P,n+l.a) 

n, A 



p,n,A J 



(ii) 



(12) 



(13) 



The heat current operator for the pure system can be determined similarly. In analogy with polarization, one defines 
the energy position operator— as 



K E = i J dr (riT(r) + H(r)r) , 



and the total Hamiltonian is H = J drH(r). Using this, one deduces the energy current from 



3 K = aR 1 



This leads to 

jE _ V F& \- 

Jx - 2 ^ 

p,a 

t e _ iv Fe \ - 
J y ~ 9 



E(Q, 



V2 
E{0,a) 

~7T 



-(c+c p ,o,a + c+ q c p ) + ^ ^{E{n + 1, a) + E(n, A))(c+„ +lct c p ,„, A + c+ n ^c Ptn+ho! ) 



n , A 



(c+c p ,o,a - c+ 0iQ c p ) + ^ -(-B(n + 1, a) + E(n, A))(c+ nA c p ,„+i :Q - c+„ +liQ c p ,„,; 



(14) 
(15) 

(16) 
(17) 
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These follow naturally from the electric current operator, after multiplying each term with the corresponding mode 
energy. Note, that the energy of the state labeled solely by (p) is zero, it belongs to the state situated at the meeting 
point of the two cones. Finally the heat current operator is related to the energy current by the simple formula: 
jQ = J E — ^tJ ? n is the chemical potential. So far we have considered the particle-hole symmetric case with /j, = 0, 
but we can easily use a finite chemical potential to break this symmetry, and introduce finite Seebeck coefficient. 



III. IMPURITY SCATTERING IN THE PRESENCE OF MAGNETIC FIELD 



In the presence of impurities, an extra term is added to the Hamiltonian: 



ff irap = KV{(r- ri ). (18) 



As a result, the explicit form of the previous operators might change. However, using Eq. (|18[) . the electric current 
remains unchanged, but the heat current changes due to the non-commutativity of the impurity Hamiltonian and the 
energy position operatoi— . As a result, impurities need to be taken into account not only in the calculation of the 
self-energy, but also in the form of the operators, and one has to use the same level of approximation for both. 

However, to avoid this difficulty, one can replace the energy terms in J E by the Matsubara frequency^, since from 
the poles of the Green's function, this will pick the appropriate energy. This replacement works perfectly in the case 
of impurities as well, when quasiparticle excitations possess finite lifetime. 

Since graphcnc is two-dimensional, positional long range order (i.e. lattice formation) is impossible at finite tem- 
peratures, since thermal fluctuations will destroy it. This is why the introduction of defect is natural in this system. 
To mimic disorder, we have chosen to spread vacancies in the honeycomb lattice, which can be modeled by taking the 
impurity strength (V) to infinity. 

To take scattering into account, we have to determine the explicit form of H imp in the Landau basis. Then the 
standard non-crossing approximation can be used 2 ^, which, in the case of graphene, is called the full self-consistent 
Born approximation due to the neglect of crossing diagram a 25 ! 26 ' 27 . Averaging over impurity positions is performed 
in the standard way. By letting the impurity strength V — > oo, which would correspond to the unitary scattering 
limit in unconventional superconductors^, we arrive to the following set of equations 2 ^: 

G(iu n ,k,n,a) = — — — — — r, (19) 

iu> n — h(n,a) — hi(ioj n ) 

G{iw n ,k) = - ^7-^, (20) 

lUJ n - L 2 (lW„) 

where 

Ei(iw„) = t-—. J^' -r, (21) 

g c [G(iuj ni k) + }_^ n a G(iw n ,k, n, a)\ 

S 2 («w„) = 7— — . \ ^ m- 1 IV ( 22 ) 

where g c = 1/(N + 1) is the degeneracy of a Landau level per unit cell, rii is the impurity concentration. The 
summation over Landau levels can be performed to yield 



G(iu n , k, n, a) = 2— {*(1 - z 2 ) - ^{N + 2- z 2 )} , (23) 



where z = (iuj n — Si(io; n ))/a; c , ^f(z) is the digamma function. The self-consistency equations can further be simplified, 
and after analytic continuation to real frequencies (iu> n — > cu + «0 + ), we can read off 

Ei(w) = S 2 (w)— — — — — -. rr . (24) 

3 c S 2 (w) + 2n,i(uj - S 2 (w)) 



At zero frequency, this simplifies to 

ti„ cm = v, cnw 1 - 

2m(N + i 



S a (0) = S 1 (0)(l- - ,1 —r ). (25) 
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The imaginary part of the self energy is always negative to ensure causality. This means that the last term in 
parenthesis on the right hand side must always be positive to assure the same sign of the imaginary parts of the self 
energies. This translates into 

Ui ^ —r — r. (26) 

2(N+ 1) K ' 

For each impurity concentration, there is a certain magnetic field strength (when TV = (l/2m) — 1), above which our 
approximation breaks down. For higher field, the self energy at zero frequency needs to be zero to fulfill Eq. (|25p and 
causality. This means, that at a finite impurity concentration, we still have excitations in the system with infinite 
lifetime. Further, we are going to show that this occurs not only on the zeroth Landau level, but on all Landau levels 
for field exceeding the critical one. To improve on this, crossing diagrams need to be considered, which is beyond the 
scope of the present work. Hence we restrict our investigation to fields allowed by Eq. (|26j) . The larger the impurity 
concentration, the larger the magnetic field we can take into account. 

The quasiparticlc density of states can be evaluated from the knowledge of the Green's function, and it reads as 

P {u) = ^Im-i-v (27) 

Without impurities, the density of states consists of Dirac-dclta peaks located at zero frequency and at E(n,a). 
By introducing impurities in the system, we expect the broadening and shift of these levels, and it can be determined 
from the solution of the self consistency equations. 

For large magnetic fields (small N), we can still solve the self-consistency equations Eqs. (|2"Tj) - ([2"2"j) . but we discover 
Dirac-delta peaks at the position of the levels and small islands between them (Fig. ^ N ~ 100). This signals that 
the non-crossing approximation is insufficient to provide these peaks with a finite broadening. As we decrease the 
field (increase N), the peaks and islands merge, and all excitations possess finite lifetime, but clean gaps are still 
observable between the levels. By further decreasing the field, the gaps disappear, the density of states becomes 
finite for all energies, and small successive bumps remain present due to Landau level formation, which tend to be 
smoothened by further decreasing the field. In this limit, the resulting density of states is very close to that in a 
d-wave superconductor—, stemming from its linear frequency dependence in the pure case. 

The broadening of the levels is not symmetric, more spectral weight is transferred to the lower energy part, which 
arises from the important energy dependence of the imaginary part of the self energies. Also the level position is 
modified in the presence of impurities due to the finite real part of the self energies, and this shift increases with the 
impurity concentration. 

The numerical solution of Eqs. (|2"Tj) - (|2"2")) . and the resulting density of states is shown in Fig. [51 From this, one can 
conjecture, that a given rij and N can qualitatively well describe different fields and concentrations, if their product 
(riiN) is the same. These features, including the non-Lorentzian broadening of the Landau levels and the development 
of small islands between the levels should be observable experimentally by scanning tunneling microscopy, for example. 

IV. ELECTRIC AND THERMAL CONDUCTIVITIES 

Using the spectral representation of the Green's functions, we can evaluate the corresponding conductivities after 
straightforward but lengthy calculations. These are related to the time ordered products of the form^i 

P 

rABr- A _ / j — iur/T jA'\ T B 



Uff(iuj) = - / dre™ (T T Jf (r)Jf (0)), (28) 



o 



where A and B denote the electric or heat current, i and j stand for the spatial component. These can be expressed 
with the use of the following transport integrals^: 

oo 



4T CO sh 2 ((e- A i)/2T) V T 
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FIG. 2: (Color online) The density of states is shown in the left panel for m = 0.001 for N = 100 (red), 1000 (blue), 3000 
(black), 10000 (magenta). The vertical red lines stand for the Dirac-delta peaks for N = 100. The right panel visualizes the 
rii = 0.01 case for N = 100 (red), 200 (blue), 300 (black). The clean case without magnetic field (N = oo) is also plotted for 
comparison in both panels (blue dashed line). 



where 



IE 



ImS 2 (e) 



ImEi(e) 



(x - ReS 2 (e)) 2 + (Im£ 2 (e)) 2 (x - E(0, a) - ReSi(e)) 2 + (ImSi(e)) 2 
ImSx(e) ImSi(e) 



+ 



2^{x- E{n,a) - Re£i(e)) 2 + (ImEi(e)) 2 (x - E(n + 1, A) - ReEi(e)) 2 + (Iml]i(e)) 2 



(30) 



is the dimcnsionlcss conductivity kernel. With the use of these, we obtain the various transport coefficients as usual: 

(31) 
(32) 

(33) 

1 L2L0 — L\ 



2e 2 

s = l T- 
- = — (l 2 -^- 



aT 



T 2 



(34) 



Here, a is the electric conductivity, S is the Seebcck coefficient, k is the heat conductivity, where the last term ensures 
that the energy current is evaluated under the condition of vanishing electric current, and L is the Lorentz number. 
Off diagonal components of the conductivity tensors, such as the Nernst coefficient, are also of prime interest, but 
they cannot be simply evaluated from Kubo formula. Even in the case of a normal metal with parabolic dispersion, 
the Kubo formula turned out to be invali d 35 ' 36 , and additional corrections have been worked out. Their determination 
for two-dimensional Dirac fcrmions is beyond the scope of the present investigation. 

For the particle-hole symmetric case {[i = 0), the Seebeck coefficient is trivially zero. If we consider the zero 
temperature, half-filled case, and assume small magnetic fields, we obtain the universal conductivity given by 



Co = 



2f_ 



and similarly for the thermal conductivity as 



K 

T 



2e 2 n 
3h ■ 



(35) 



(36) 
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The Seebeck coefficient is zero. From this, the Lorentz number takes its universal value 



(37) 



which means, that in this limit, the Wiedemann-Franz law hold a 16 i 23 . Landau levels always develop around the 
meeting point of the conical valence and conduction band. If we are at half filling (/i = 0), no levels cross « when 
varying the magnetic field, since they are symmetrically placed below and above. However, when u is finite, Landau 
levels can cross its value with changing the field, and we expect Shubnikov-dc Haas oscillations. In general, when 
the number of levels below /i is large (or oj c <C |/z|), we can conjecture the periodicity of these oscillations. Assume 
that a level (the nth) sits right at the chemical potential (fi = uj c \Jn + 1). Then, the distance from the adjacent level 
determines the period of the oscillations. This is 

\w( ±1 ^ w( m u c w c v 2 F e\B cos(9)\ 

\E(n + l,a)-E(n,a)\K —== = — = B (38) 

2\Jn + 1 2(i fx 

provided, that n^$> 1. This means, that albeit the Landau levels show an unusual oc ^fri dependence of the level index 
compared to that in a normal metal oc n, the Shubnikov-de Haas oscillations turn our to be still periodic as a function 
oil/ B. The comparison of the coefficient of the magnetic field in Eq. (|3"5)l to that in a parabolic band2& suggests, that 
the cyclotron mass can be defined as m c = ii/v F . Even though the spectrum is linear, the finite chemical potential 
provides us with an energy scale for mj±. This can readily be checked in Fig. [31 where not only the field, but the 
angle dependence of the conductivity of shown for different field strength. The larger the magnetic field, the more 
visible the oscillations are, although these can be smeared by increasing the concentrations. 

The explicit value of the chemical potential, which is fixed by the particle number at a given temperature and 
field, should also be determined self-consistently. However, no serious deviations from its initial values have been 
detected during the evaluation process, and these did not affect the dependence of physical quantities on T and B 
in the investigated range of parameters. Presumably, taking a large value of the chemical potential would require its 
self-consistent determination as well. 

In Fig. [U we show the magnetic field dependence of the heat conductivity. It resembles closely to the electric one at 
low temperatures. However, at higher temperatures, each peak in the oscillations split into two. This occurs, because 
in the electric conductivity, the kernel is sampled by the 1/ cosh 2 ((e — fx)/2T) function, which gathers information 
about excitations at the chemical potential. However, an extra (e — /i) 2 factor appears in the heat response, which 
measures the immediate vicinity of u above and below, within a window 2T, which gives the splitting. The oscillations 
become smoothened with decreasing field, in contrast to Ref. [l8l where large oscillations were found even at small 
fields. The difference can be traced back to our field dependent scattering rate (Eqs. pTjl - P^)) ). as opposed to the 
field independent one used in Ref. Similar features have been observed in highly oriented pyrolytic graphite"*. 
By decreasing the field, N increases, and the density of states becomes similar to that of a d-wave superconductor—, 
without significant deviations from linearity. Both a and k decrease with field, a feature already present at fx = 0. As 
we increase the field, oj c increases, and so does the distance between Landau levels. Then, at a given temperature, a 
smaller number of states will be present for excitations around fx, hence the corresponding conductivity decreases. The 
Seebeck coefficient shows sharp oscillations which die out with temperature. Its background value, after subtracting the 
oscillations, is found to be almost magnetic field independent, but smoothly increases with temperature. The Lorentz 
number remains close to one, if we subtract the oscillations. However, due to the double (single) peak structures 
in the heat (electric) response, their ratio shows wild, but sharp deviations from unity at specific fields, where the 
Wiedemann-Franz law is violated. In contrast to this, one would have encountered large and wide oscillations in the 
Lorentz number as a function of field in the presence of phenomenological, constant scattering rate. 

In Fig. [5j we show the evolution of the electric and heat conductivity and the Seebeck coefficient as a function 
of chemical potential. In accordance with experiment in Ref. 0, we also find oscillations, corresponding to Landau 
levels, which also smoothen with temperature. Interestingly, the splitting of the peaks in the heat conductivity is 
nicely observable as a function of fx. These occur in such a way that they produce antiphase oscillations with respect 
to the electric one, and lead to the violation of the Wiedemann-Franz law. The Seebeck coefficient shows peculiar 
behaviour. At the particle-hole symmetric case, it is zero, and remains mainly so apart from large oscillations. 

The temperature dependence of the electric and heat response is shown in Fig. [SJ Both increase steadily with 
temperature, since more available states are accessible with T. However, at small temperatures, a small decrease is 
observable in low fields, in accordance with other studie a 16 i 25 The Seebeck coefficient first increases, and after a broad 
bump, decreases with T. For higher temperatures, the bandwidth D makes its presence felt. The Wiedemann-Franz 
law remains intact at low temperatures and fields, but becomes violated for higher T or B. 
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FIG. 3: (Color online) The angular dependent magnetoconductivity oscillations are visualized for fi — 0.05D, 7ii=0.001 and 
T = 0.0001D, for magnetic fields N = 600 (red), 1000 (blue), (2000) (black), (3000) (green) and 5000 (magenta) in the left 
panel. With increasing field (decreasing N), the oscillations become more pronounced, signaling the discrete Landau level 
structure. The right panel shows the electric conductivity for /i = 0.05D, ni=0.001 and T/D = 0.0001 (red), 0.001 (black) and 
0.01 (blue). For higher field, we arrive to the region, where crossing diagrams need to be taken into account. The inset shows 
the electric conductivity as a function of l/|Bcos(0)| to emphasize its periodicity. 
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FIG. 4: (Color online) The left panel shows the Lorentz number as a function of the inverse magnetic field to sress the periodic 
violation of the Wiedemann-Franz law for fj, = 0.05D, ni=0.001 and T/D = 0.0001 (red), 0.001 (black) and 0.01 (blue). The 
right panel shows the heat conductivity and the Seebeck coefficient (inset) for the set of same parameters. 



V. CONCLUSION 



We have studied the effect of localized impurities in two dimensional Dirac fermions in the presence of quantizing, 
arbitrarily oriented magnetic field. The energy spectrum depends on the level index as oc y/n, as opposed to the 
n + 1/2 linear dependence in normal metals^. Expressions for both the electric and heat current in the presence 
of magnetic field were worked out. The self-energy in the full Born- approximation obeys self-consistency conditions, 
resulting in important magnetic field and frequency dependence of scattering rate and level shift. In the density of 
states, only a small island shows up close to zero frequency for small fields, similarly to d-wave superconductors^. 
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FIG. 5: The electric (blue solid line) and heat (red dashed line) conductivity and the Seebeck coefficient (inset) are shown 
as a function of the chemical potential for T = 0.001-D, TV = 1000, n; = 0.001. Due to the antiphase oscillations, the 
Wiedemann-Franz law is violated. 




T/D T/D 

FIG. 6: (Color online) The electric (blue solid line) and heat (red dashed line) conductivity are shown in the right panel for 
rii=0.001, (j, = 0.05D, iV=600, 1000, 3000 and 10000 from bottom to top. Note the 1/tt 2 reduction of the heat conductivity. 
The right panel shows the Seebeck coefficient and the Lorentz number (inset) for the same parameters from top to bottom, 
with dashed line for iV = 600. Note the violation of the Wiedemann-Franz law at low temperatures at high fields (smaller AQ! 

By increasing the field, oscillations become visible, corresponding to Landau levels. By further increasing the field, 
these become separated from each other, and clean gaps appear between the levels, in which intragap states, small 
islands show up at high field. The non-Lorcntzian broadening of Landau levels and the intragap features differ from 
previous studies assuming a constant scattering rate, and should be detected experimentally in graphenc. 

Both the electric and thermal conductance shows Shubnikov-dc Haas oscillation in magnetic field, which disappear 
for small fields and higher temperatures. These are periodic in l/B, similarly to normal metals, in spite of the 
different Landau quantization. The Seebeck coefficient shares these features, but its oscillations are really large as 
opposed to a and re. The Wiedemann-Franz law stays close to unity, except at certain fields, where large deviations 
are encountered, which vanish with decreasing field. Besides oscillations, both a and re decreases with field, since the 
larger the cyclotron frequency, the smaller the probability of finding states around fi. These are in agreement with 
experiments on the thermal conductivity of highly oriented pyrolytic graphite^. Oscillations are also present as a 
function of chemical potential, similarly to experimental findings 3 . 
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The temperature dependence of the conductivities is rather conventional, both a and k increases with temperature 
steadily, regardless to the value of the chemical potential. The Scebcck coefficient exhibits a broad bump around 
T ~ fx, and decreases afterwards. The Wicdcmann-Franz law is obeyed for small temperatures and field, but violated 
for higher values. 
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